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Abstract. Massive stars lose a large fraction of their original mass over the course of their 
evolution. These stellar winds shape the surrounding medium according to parameters that 
are the result of the characteristics of the stars, varying over time as the stars evolve, leading 
to both permanent and temporary features that can be used to constrain the evolution of the 
progenitor star. 

Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the 
characteristics of the circumstellar medium (CSM) should be observable in the signal of 
GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the 
medium it collides with, and indirectly because the GRB can only be observed through the 
extended circumstellar bubble that surrounds each massive star. 

We use computer simulations to describe the circumstellar features that can be found in 
the vicinity of massive stars and discuss if, and how, they may appear in GRB observa- 
tions. Specifically, we make hydrodynamical models of the circumstellar environment of a 
rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor 
candidate model. 

The simulations show that the star creates a large scale bubble of shocked wind material, 
which sweeps up the interstellar medium in an expanding shell. Within this bubble, tempo- 
rary circumstellar shells, clumps and voids are created as a result of changes in the stellar 
wind. Most of these temporary features have disappeared by the time the star reaches the 
end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star 
in a high density environment simplifies the evolution of the CSM as the more confined 
bubble prohibits the formation of some of the temporary structures. 
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1 . Introduction (> 12 M ) these winds can be very strong due 

to their high luminosity, which drives the wind. 
As stars evolve, they lose mass in the form As a stellar wind expands into the surround- 
of stellar wind. In the case of massive stars ing interstellar medium (ISM) it creates a large 

(> 10 pc) bubble, with a complicated internal 
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morphology reflecting the characteristics of the 
stellar wind, which in turn reflect the evolu- 
tion of the progenitor star. The morphology of 
the circumstellar bubble can be observed di- 
rectly (e.g. in the form of circumstellar nebula), 
but can also show up in observations of long 
Gamma-Ray Bursts (GRBs), which can occur 
once the central star reaches the end of its life. 
This happens in two ways: 

1 . In a GRB afterglow circumstellar gas may 
be observed as blue-shifted absorption fea- 
tures in the spectrum of the GRB afterglow 
(e.g. Prochaska et al. 2008). 

2. The circumstellar gas may influence the 
evolution of the GRB itself. As the GRB 
jet expands, it will collide with circumstel- 
lar gas, which, depending on the density of 
the circumstellar medium (CSM) can cause 
a change in the expansion velocity and lead 
to a change in the observed signal as kinetic 
energy is transferred to heat and then radi- 
ated away (e.g. van Eerten et al. 2010). 

In this paper we describe some of the features 
found in the circumstellar medium of massive 
stars and discuss if, and how, they may show 
up in the observations of GRBs. 



2. Morphology of circumstellar 
bubbles 

2.1. General shape 

The general (semi-permanent) shape of wind- 
blown bubbles was described analytically by 
Weaver et al. (1977). Close to the star the wind 
expands in all directions, leading to a 1/r 2 
density profile. Eventually this wind passes 
through the wind termination shock and en- 
ters the "shocked wind bubble", which has a 
(nearly) constant density. Since the kinetic en- 
ergy of the wind is turned into heat at the ter- 
mination shock, the shocked wind has a very 
high (~ 10 7 K) temperature, which gives it a 
high thermal pressure, causing it to expand out- 
wards, sweeping up a shell of shocked inter- 
stellar gas. 



2.2. Temporary features 

As the star evolves, changes in the massloss 
rate and wind velocity cause shells to appear 
inside the bubble. These temporary features, 
which can be observed as circumstellar nebu- 
lae, have only a relatively short (~ 10 4 — 10 5 yr) 
live span, before they dissipate. Therefore, they 
can be used to determine the evolution of the 
progenitor star and, if observed, can be power- 
ful indicators as to what type of star produces 
a certain type of supernova or GRB. 

Due to their more complicated morphol- 
ogy the evolution of such shells can best be 
described through numerical simulations. Such 
simulations have been made for the transi- 
tion from red supergiant (RSG) to Wolf-Rayet 
(WR) star (Garcfa-Segura et al. 1996; van 
Marie et al. 2005; Freyer et al. 2006) and lumi- 
nous blue variable (LBV) to WR star (Garcfa- 
Segura et al. 1996; Freyer et al. 2003; van 
Marie et al. 2007). Recent work by Toala & 
Arthur (2011) described similar interactions 
and included additional physics in the form 
of radiative transfer and thermal conduction, 
while Dwarkadas (2007) showed the interac- 
tion between a supernova and the wind bub- 
ble in which it originates. The effect of an 
anisotropic stellar wind, resulting from stellar 
rotation was included by Chita et al. (2008) and 
van Marie et al. (2008). 

3. Numerical models 

3.1. Stellar evolution models 

As input models for our numerical simula- 
tions we use one of the rapidly rotating, chemi- 
cally near-homogeneous stellar evolution mod- 
els from Yoon et al. (2006). These models fol- 
low special evolutionary tracks, which differ 
from the classical evolution due to their rapid 
rotation. Rotationally induced mixing allows 
the star to maintain hydrogen core burning un- 
til all hydrogen in the star is exhausted before 
it makes the transition to helium burning. As 
a result the star has no giant phase, but makes 
a direct transition from main-sequence to WR 
phase. This transition (which occurs at the end 
of hydrogen burning) is marked by a short pe- 
riod during which the star rotates at critical ve- 
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Fig. 1. Hertzsprung-Russel diagram of a rapidly ro- 
tating 16M Q star. The star has no giant stage and 
stays at the high temperature end of the diagram. 
The start and end of the first rapid rotation phase are 
marked, as is the start of the second rapid rotation 
phase. 
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Fig. 2. Mass loss rate and wind velocity (averaged 
over azimuthal angle) for the stellar evolution model 
from Fig. 1. At the pole, massloss rate and wind ve- 
locity are higher than the average. At the pole, they 
are lower. 



locity. Because such stars maintain a high ro- 
tation rate until the end of their existence, they 
are thought to be GRB progenitors. One such 
model was used in van Marie et al. (2008). We 
now use the A16f0.5 model from the same pa- 
per to explore in more detail the interactions 



that take place in the circumstellar environ- 
ment. This model represents a 16 M star with 
an initial metallicity of 0.002 and an initial ro- 
tation rate of 50% of critical rotation. The evo- 
lutionary track of this star is shown in Fig. 1. 
From this stellar evolution model we obtain 
massloss rate, wind velocity and rotation rate 
as a function of time (Fig. 2), which are used as 
input for the hydrodynamical simulations. The 
episode of critical rotation, which marks the 
end of core hydrogen burning, is clearly visible 
as a massive increase in massloss rate at about 
1.812 x 10 7 yrs. Afterward, the massloss rate 
decreases again. A second outburst starts at the 
end of the stellar evolution (~ 1.86 x 10 7 yrs), 
but does not reach a similar height. During the 
outbursts, wind velocity is severely reduced. 

3.2. Numerical hydrodynamics 

For our hydrodynamical simulations we use 
the MPI-AMRVAC code (Keppens et al. 2011). 
This code solves the conservation equations of 
mass, momentum and energy on an adaptive 
mesh grid and allows us to choose from a wide 
variety of solvers. We include optically thin ra- 
diative cooling, using the table from Schure 
et al. (2009), which was generated with the 
SPEX package (Kaastra et al. 1996). 

For the main sequence, during which no 
major changes in stellar wind parameters take 
place, we run the simulation in 1-D to safe 
computation time. We use a spherically sym- 
metric grid with a maximum radius of 140 pc 
and a basic grid of 400 points. We allow 6 
additional levels of refinement, so that our 
grid effectively has 25 600 points. This grid is 
filled with a constant density of either 2 par- 
ticles per cirT 3 (Simulation A), representative 
of the general interstellar medium, or 2 000 
particles per cirT 3 (Simulation B), represen- 
tative of the dense molecular clouds in which 
massive stars form. At the inner radial bound- 
ary the grid is filled with matter according to 
the time-dependent stellar wind parameters ob- 
tained from the stellar evolution model. 

As the star approaches the critical rotation 
phase, it becomes necessary to model the cir- 
cumstellar bubble in 2-D. Therefore, we map 
the 1-D result onto a 2-D grid. For this we 
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Fig. 3. Density of the circumstellar medium at the 
end of the main sequence (after 1.81 x 10 7 years) 
for both simulations A and B, just before the star 
reaches critical rotation. The solid line represents 
the simulation with high density ISM. 



choose a spherical grid in the r-8 plane, with 
a radius of 97.5 pc and a 90° opening angle 
starting at the polar axis. The basic grid has 
300x200 points and is allowed two more lev- 
els of refinement, giving us an effective grid of 
1200x800 points. The azimuthal dependence 
of the massloss rate and wind velocity are 
calculated according to the gravity darkening 
model by Dwarkadas & Owocki (2002), 

M oc 1-Q 2 sin 2 6» (1) 
sin 8 (2) 

with M the massloss rate v w j n d the radial wind 
velocity,Q = v/v cr j t the rotation rate, relative to 
break-up speed and 6 the co-azimuthal angle. 
This model provides a more accurate descrip- 
tion for radiatively driven winds than the older 
"wind compressed disk" model (Bjorkman & 
Cassinelli 1993) used by Chita et al. (2008) 
and van Marie et al. (2008). Because the grav- 
ity darkening model is only valid for the stellar 
wind and does not take into account the effect 
of matter being driven of the equator by cen- 
trifugal forces we truncate the rotation rate to 
99% of critical rotation. 
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Fig. 4. Circumstellar density for simulation A ~ 
60 000 years after the result in Fig. 3, just after the 
critical rotation phase. During the critical rotation 
phase, which channels most of the wind toward the 
pole, the high massloss rate and low wind velocity 
cause the formation of a bipolar shell. Afterward, he 
wind velocity increases, leading to the formation of 
a third shell. 
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Fig. 5. Similar to Fig. 4 at the end of the stellar evo- 
lution. The shells have disappeared, leaving a highly 
turbulent wind bubble. 

4. Results 

4.1. 1-D models 

At the end of the 1-D simulation (1.81 X 
10 7 years after the start of the evolution), just 
before the star reaches critical rotation, the 
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CSM morphology conforms to the analyti- 
cal description from Weaver et al. (1977), as 
shown in Fig. 3. The effect of a high density 
ISM is clearly visible, as the bubble of simu- 
lation B is far more compressed. This effect, 
which is predicted by the analytical approxi- 
mation and shown numerically in van Marie 
et al. (2006), is enhanced by radiative cooling, 
which is far more effective at high densities, 
causing energy to be lost from the system, the 
radiative cooling also causes the internal struc- 
ture of the swept-up shell, which has the high- 
est density at the contact discontinuity. As mat- 
ter is swept up into the shell it starts to cool. 
The further it travels into the shell, the more it 
cools, causing to to reach its lowest tempera- 
ture (and therefore highest density) at the con- 
tact discontinuity. 

4.2. 2-D models 

As the star in Simulation A reaches critical 
rotation, the wind velocity decreases and the 
massloss rate increases causing the formation 
of a bipolar shell, which moves outward into 
the hot bubble. The shape of the shell results 
from the azimuthal dependence of the wind, 
which directs most of the mass towards the 
pole (Dwarkadas & Owocki 2002). Once the 
critical rotation phase ends, the wind veloc- 
ity increases once again, causing the wind to 
sweep up its surroundings in a new shell, which 
will eventually overtake its bipolar predecessor 
(Fig. 4). This sequence of events in the CSM 
shows a close resemblance to the results from 
Chita et al. (2008) and van Marie et al. (2008), 
despite the use of a completely different de- 
scription for the wind of a rotating star. At the 
end of the stellar evolution (Fig. 5), both shells 
have disappeared after colliding with the outer 
shell and dissipating into the hot bubble. They 
leave behind a very turbulent region, especially 
over the pole of the star. 

Figures 6-7 show the CSM for simula- 
tion B 40 000 years after the end of the 1-D 
simulation and at the end of evolution respec- 
tively. Because of the much smaller bubble, the 
bipolar shell that forms during critical rotation 
collides very early with the outer shell and dis- 
appears. The third shell never even has time to 
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Fig. 6. Circumstellar density for simula- 
tion B ~ 40 000 years after the result in Fig. 3. 
Because of the small size of the bubble, the bipolar 
shell has already collided with the outer shell in the 
polar region. At the equator the bipolar shell is still 
moving outward. The third shell never forms. 




Fig. 7. Similar to Fig. 4 at the end of the stellar evo- 
lution. Unlike in simulation A the hot bubble shows 
very little turbulence. 



form. As a result, the CSM is much less turbu- 
lent at the end of the star's life. 

Although the star spins up again at the end 
of its evolution, it does not again reach critical 
rotation and this episode has only minimal in- 
fluence on the CSM, unlike the 20 M model 
described in van Marie et al. (2008) 
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5. Discussion and Conclusions 

The density of the ISM plays an important role 
in the evolution of the circumstellar bubble. 
The high density medium constrains the radius 
of the bubble, which in turn limits the forma- 
tion of temporary features, such as circumstel- 
lar shells. Th effect of these limitations still 
shows at the end of the star's life, when the 
high density simulation shows much less tur- 
bulence than the more extended bubble in the 
low density medium. 

The temporary shells in the CSM, created 
by the changes in the stellar wind during the 
critical rotation phase, have disappeared by the 
time the star reaches the end of its evolution 
and the second phase was not powerful enough 
to generate a similar event. The critical rota- 
tion phase has left behind a great deal of tur- 
bulence in the shocked wind bubble. This re- 
sults in density differences of more than order 
of magnitude that may influence the emission 
signal of the GRB afterglow as it sweeps up the 
high density regions. 

The absorption spectrum generated by 
this bubble should be unremarkable as the 
shells have disappeared, leaving only the free- 
streaming wind to generate a blue-shifted ab- 
sorption feature (van Marie et al. 2005, 2007, 
2008) and this is so close to the star (and 
therefore the GRB), that only very high ioniza- 
tion states can survive (Prochaska et al. 2007, 
2008). The outer shell might be visible as an 
independent absorption feature, but its veloc- 
ity is low (~ 25 km/s), making it nearly indis- 
tinguishable from the ISM. 
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